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ABSTRACT 

Context. Central galaxies in galaxy clusters may be key discriminants in the competition between the cold dark matter (CDM) 
paradigm and modified Newtonian dynamics (MOND). 

Aims. We investigate the dark halo of NGC 1399, the central galaxy of the Fornax cluster, out to a galactocentric distance of 80 kpc. 
Methods. The data base consists of 656 radial velocities of globular clusters obtained with MXU/VUT and GMOS/Gemini, which is 
the largest sample so far for any galaxy. We performed a Jeans analysis for a non-rotating isotropic model. 

Results. An NFW halo with the parameters r, = 50 kpc and p, = 0.0065 Mo/pc^^ provides a good description of our data, fitting 
well to the X-ray mass. More massive halos are also permitted that agree with the mass of the Fornax cluster as derived from galaxy 
velocities. We compare this halo with the expected MOND models under isotropy and find that additional dark matter on the order of 
the stellar mass is needed to get agreement. A fully radial infinite globular cluster system would be needed to change this conclusion. 
Conclusions. Regarding CDM, we cannot draw firm conclusions. To really constrain a cluster wide halo, more data covering a larger 
radius are necessary. The MOND result appears as a small-scale variant of the finding that MOND in galaxy clusters still needs dark 
matter. 
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1. Introduction 

The dark matter halos of early-type galaxies are investigated 
less well than those of spiral galaxies (for a recent review, see 
Romanowsky 2006 ). Globular clusters (GCs) and planetary neb- 
ulae are the main stellar dynamical tracers, while X-rays can 
be used at even larger galactocentric distances for a gas dy- 
namical approach. The reliability of a deduced mass profile de- 
pends strongly on the number of dynamical probes available. 
Regarding GCs, it is therefore natural that the rich GC systems 
of bright elliptical galaxies have so far been the preferred targets 
for such stu dies, f or example M87 (Cote et al.l ^UOTT l , NGC 4472 
(Cote et al. l2003T l. and NGC 4636 (Schuberth et al. l2006] l. The 
largest sample until today has been measured for NGC 1399, the 
central galaxy of the Fornax galaxy cluster. Richtler et al. (120041 
Paper I) analyzed the velocities of about 450 GCs out to a galac- 
tocentric distance of 40 kpc and found an approximately con- 
stant velocity dispersion. The coiTesponding dark matter halo 
was found to be consistent with an NFW-halo (see section 5 for 
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the numerical values). The relatively poor knowledge that we 
have until now regarding the dark halos of early-type galaxies, 
does not indicate a common behavior such as the constant ro- 
tation curves found in spiral galaxies. While the inferred cir- 
cular velocities of NGC 1399 or NGC 4636 appear to be flat, 
they seem to rise for M87, and most strikingly so for NGC 6166 
(Kelson et al.'2002^, the central galaxy in Abell 2199. In ellip- 
tical galaxies of lower luminosity, the projected velocity disper- 
sion can even be falling (Romanowsky et al. I2003I I. perhaps in- 
cluding the circular velocity. One may speculate that these relate 
to the environment. The dark halos of central galaxies might be 
determined by the dark halos of their host galaxy clusters rather 
than being halos belonging to individual galaxies. However, a 
larger sample of well-studied galaxies is necessary to arrive at 
firm conclusions. 

We have now augmented the number of GC velocities around 
NGC 1399 to nearly 660 reaching a galactocentric distance of 
80 kpc, which is an unprecedently large and extended sample. 
In this Letter, we present the implications of this sample for 
the dark halo, updating the conclusions of Paper! Moreover, 
we discuss the dark halo in the context of Modified Newtonian 
Dynamics (MOND) (e.g. Milgrom 1983 Sanders & McGaugh 
2002) . The detailed description of the observations, the reduc- 
tion, and the analysis will be presented in a forthcoming paper 
(Schuberth et al., in prep), taking into account the population 
substructure of the GC system as well as more anisotropic mod- 
els. 
In the following, we adopt a distance of 19 Mpc for NGC 1399. 
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2. Observations 

Our new data were obtained with the multi-object- 
spectrographs FORS2/MXU and GMOS at the VLT and 
Gemini-South observatories, respectively. The GC candidates 
were selected from the wide-field photometry by Dirsch 
et al. (I2003I I. The spectroscopic observations of ten MXU 
masks in seven fields were carried out in visitor mode during 
three nights (1-3 December 2002) at the European Southern 
Observatory Very Large Telescope (VLT) facility on Cerro 
Paranal, Chile. The observations with the Gemini Multi-Object 
Spectrograph (GMOS) on Gemini-South were carried out in 
queue mode in November 2003 and December 2004. A total of 
ten spectroscopic masks in five fields were observed. We then 
merged the velocity catalogues with the previous velocities of 
Paper I after having revised the old spectra by cross correlation 
methods with the same templates used for the new data. The 
final velocity database comprises 656 GCs; 160 new targets 
come from FORS2 and 73 targets from GMOS-S. A typical 
velocity uncertainty is 50 km/s, while it can be considerably 
higher for the faintest targets, which have a brightness of about 
niR X 22.8. 

3. Velocity dispersions 

Converting measured velocities to velocity dispersion is not 
straightforward. Fluctuations may arise from inhomogeneous 
sampling or from small numbers of outliers (stemming e.g. 
from erroneous velocities or contaminants belonging to other 
galaxies — in our case NGC 1404, which is at a projected dis- 
tance of 54 kpc). Therefore we define a "safe" subsample, 
where the objects have distances greater than 3 arcmin to 
NGC 1404, velocity uncertainties smaller than 75 km/s, and 
magnitudes in the range 21 < m^ < 22.35. We calculate disper- 
sions by applying the maximum likelihood estimator of Pryor 
& MeylanU993), with the individual velocity uncertainties as 
weights. For the full sample, the dispersion was calculated by 
moving bins, each bin containing 60 objects with a step of one 
object. For the "safe" sample, we calculate the dispersion in in- 
dependent bins. The sample velocities and derived dispersions 
are shown in Fig.[T] 

The velocity dispersion is approximately constant out to 
80 kpc, thus continuing the trend found in Paper I. The various 
curves refer to different models, which are explained in the fol- 
lowing sections. 

4. l\/lodels 

In Jeans models that only use velocity dispersions and not higher 
moments of the velocity distribution, the mass distribution and 
the orbital anisotropy of the tracer population are degenerate. 
Here we mainly consider isotropic models, which are good ap- 
proximations at least to the entirety of those clusters that were 
found in Paper I to best explain the GC dispersions (see also Cote 
et al. 1200 11 1. In the MOND discussion, we also give fully ra- 
dial models (which are not realistic) as a guide to how a radial 
anisotropy affects the results. 

We model the projected velocity dispersions under isotropy 
according to 
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where N(R) is the projected number density of the tracer popula- 
tion, (Tios the line-of-sight velocity dispersion to be compared to 
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Fig. 1. Kinematics of GCs around NGC 1399. The upper panel 
shows the radial velocities of our total sample versus the ra- 
dial distance, where 80 kpc correspond to 14.5 arcmin. The two 
crosses are outlier velocities that we have discarded, and the 
larger points show objects from the "safe" sample. The horizon- 
tal line marks the systemic velocity, and the vertical dashed line 
indicates the radial limit used for the "safe" sample. The middle 
panel shows the velocity dispersion profile. The thin solid line 
with shaded region shows the result derived for the full GC sam- 
ple along with its uncertainties, while the circles with error bars 
show the "safe" sample. The curves represent various dynamical 
model fits (isotropic), as described in the text, and with param- 
eters summarized in the legend. The lower solid line represents 
the expectation from MOND under isotropy. The lower panel 
shows in addition fully radial MOND models (without gas) with 
a cut-off radius of 250 kpc. 



our observed values, r is the radial distance from the center and R 
the projected distance, £ the spatial (i.e., three-dimensional) den- 
sity of the GCs, and M{r) the enclosed mass. See e.g., Mamon 
& Lokas (2005^ and van der Marel ( 119931 ) for more general ana- 
lytical solutions of the Jeans equation. 

The surface number densities of GCs are adopted from 
Bassino et al. (2006) as the sum of the number densities for 
metal-rich and metal-poor clusters, resulting in the 3D density 
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with the values N^ = 35.54, /^o = lf74, a = 0.84, and the Beta- 
function S (i, a) = 2.22. 

As a model for the luminous matter, we adopt the power 
law given by Dirsch et al. (2003) resulting from the wide-field 
photometry of NGC 1399 together with a (constant) stellar mass 
to light ratio, which we fix to MjL^ - 5.5, following Paper I. 
This is about the maximum value, which is found for early-type 
galaxies (see also SagHa et al. l20U0l Bell et al. 12003]) . 

For the shape of the dark halo, we assume an NFW type 
halo, which for the inner 40 kpc already gave a good description 
(Paper I). The cumulative mass is 
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where g^ and r^ are the characteristic density and radius, respec- 
tively. 



5. Mass profiles 

The full sample can be reproduced well by adding an NFW halo 
with Qs - O.OO85M0/pc^ and r, - 50 kpc, the "safe" sample 
with Qsfl - O.OO65M0/pc^ and r, = 50 kpc (long-dashed and 
dashed fine in Fig. [B to the mass of NGC 1399. As TableQ 
shows, the gas mass (which is included in the halo) is negligible 
within 80 kpc with respect to the stellar mass and unimportant 
within the virial radius with respect to the dark halo. 

The uncertainties are mostly of a systematic nature and dif- 
ficult to quantify, but these two halos probably embrace the ac- 
ceptable values with a preference near the "safe" sample in the 
sense of the best-fit halos. 

These halos can now be compared with the X-ray mass pro- 
file presented by Paolillo et al. (2002). Figure|2] shows the total 
mass for the various models; i.e. the NFW halo includes dark 
matter, gas, and the stellar matter outside NGC 1399; however, 
it is vastly dominated by dark matter. The overall agreement of 
the X-ray mass profile, particularly with the "safe" halo, is very 
good. There is no support however for the steep rise between 
20 and 40 kpc and the subsequent flat part, which would mean 
a strong decline in the mass density to almost zero. This feature 
has been already noted in the ASCA study of Ikebe et al. ( 1996), 
who interpreted it as the transition radius between the galaxy 
and the cluster potential. Such a mass profile would probably be 
dynamically unstable, so we suspect that it is possibly a conse- 
quence of azimuthally averaging over a disturbed intensity field. 
The "safe" halo has a high concentration and only marginally 
fits to the simulations of Bullock et al. (120011 ). but better than the 
halo from Paper I. The virial radius is 670 kpc, the mass inside 
this radius is 1.8 x IO'-'M©, and the concentration c = 13.4. Until 
now, the only dynamical analysis of the Fornax cluster is that of 
Drinkwater et al. (2001. hereafter DEA), who used galaxy ve- 
locities to constrain the cluster mass. 

To investigate what one would expect from the DEA mass for 
the projected velocities in our radial range, we describe the 
DEA mass by an NFW profile (which includes dark matter 
and baryons, but not NGC 1399) with the parameters QsXi - 
O.OOO5M0/pc^ and r, = 220 kpc, which reproduces the solid 
line in Fig. 4 of DEA out to 1 Mpc very well. To compare this 
model to the GC dispersions in the inner regions of the halo, we 
add in the stellar mass of NGC 1399, and show the modeling re- 
sult in Fig.[T] Although this halo model was not among our best 
fits, it is also not excluded, as it seems to follow the velocity 
dispersions quite well. However, it falls below the X-ray mass 



Table 1. The table lists the masses of different components or 
halos for two radii, the larger one being the virial radius of the 
"safe" halo. All values are given in units of IO'^Mq. The stellar 
mass at 670 kpc is an extrapolation of the NGC 1399 luminos- 
ity profile and serves as a proxy for the total stellar mass. The 
gas mass is taken from Paolillo et al. 2002 and is extrapolated 
beyond 200 kpc. The "safe" halo is explained in the text. "DEA- 
mass" denotes the Fornax cluster mass according to Drinkwater 
etal. l200n 
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Fig. 2. Comparison to X-ray measurements. The grey area 
shows the range of the mass models presented by Paolillo et 
al. (2002). The short-dashed line refers to our preferred halo, 
the long-dashed line to the halo derived from the full sample. 
The solid line is the halo derived in Paper I. The stellar mass 
is shown as dashed-dotted line. The dynamical MOND mass is 
shown as the dot-dot-dashed line. 

model at small radii (Fig. |2]). The difference between the mass 
profile of DEA and our "safe" halo is listed in Table[T] The gas 
mass of Paolillo et al. (2002) is also given. It is only known out 
to about 200 kpc, and thus is extrapolated to larger radii. Table[T] 
shows that the agreement of our halos with the mass of DEA is 
unsatisfactory. It becomes even worse at larger radii: at 1 Mpc, 
the DEA mass is 6 ± 2 x IO'^'Mq versus our mass of 2.2 x IO'^Mq 
(2.9 X IO'^Mq for the more massive halo). We conclude that, al- 
though our sample is already the largest one available, we are 
not able to conclusively distinguish between different halos on a 
cluster-wide scale. There is, however, evidence from the data of 
Bergond et al. (2007) that our halo might be extrapolated out to 
200 kpc (Schuberth et al. 12007] ). 

6. Modified Newtonian dynamics 

As an alternative to the CDM paradigm, Milgrom's MOND (e.g. 
Milgrom 119831 Sanders & McGaugh 120021 ) is mostly success- 
ful in accounting for the kinematics of disk galaxies and/or the 
baryonic Tully-Fisher relation (McGaugh 2005). Little is known 
about early-type galaxies. Schuberth et al. (2006) find that 
NGC 4636 is indeed consistent with MONDian kinematics. Also 
the declining velocity dispersion of NGC 3379 (Romanowsky 
et al. 120031) does not necessarily contradict MOND (Milgrom 
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& Sanders I2003I I. However, on the scale of galaxy clusters, 
MOND apparently needs additional dark matter (e.g. Sanders 
[200 3 Pointecouteau & Silk l2005T l. It is therefore interesting to 
investigate whether NGC 1 399, as a central galaxy, is compatible 
with MOND or not. 

Assuming spherical symmetry, we apply the classic MOND 
recipe where g and g^ are the MOND acceleration and 
Newtonian acceleration, respectively, aq a universal constant, 
and fi{g/a^)) is a function interpolating between the Newtonian 
and the MONDian regimes. 

Following Famaey & Binney ( 120051) . we adopt// = x/(l +x) 
with X - g/ao (see also Zhao & Famaey |20061 l. The circular 
velocity curve then reads 
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where WdrcN is the Newtonian circular velocity. For ao we 
adopt the value recommended by Famaey et al. (120071 ): 1.35 x 
10 '^cms"^. The analysis strongly depends on the shape of 
the interpolation, since g is comparable to ao even at 80 kpc. 
Adopting the Bekenstein ( |2004| l interpolation function would in- 
crease Vcirc significantly, but it does not fit spiral rotation curves. 

To calculate the MOND projected velocity dispersion, we 
used the standard Jeans equation, with the Newtonian mass cor- 
responding to the MONDian circular velocity. The MOND mod- 
eling results are shown in Fig.[T] Under isotropy, the MOND dis- 
persions are significantly lower than the profile derived from the 
GC dynamics. 

Fully radial models are shown in the lower panel. We caution 
that, in contrast to the isotropic case, the outer integration limit 
now gains importance. Integrating to infinity, the MOND dis- 
persion could reproduce the data. However, we do not have in- 
finitely distant GCs. The extension of the GCS is about 250 kpc 
(Bassino et al. '2006V The model for cut-ofFs of 250 kpc shows 
that a more realistic radial model does not help in achieving bet- 
ter agreement with MOND. We therefore conclude that MOND 
does not explain the dynamics of NGC 1399 without additional 
dark matter 

Such an additional hypothetical dark halo must have a rela- 
tively small core and a rapidly declining density profile in order 
to avoid an increasing circular velocity. Without claiming pre- 
cision, a halo like Qadd = 0.05Mo/pc''(l + r/12)""^, where r is 
in kpc, would suffice. It corresponds to adding a dark mass of 
7 X 1O"M0 within a radius of 80 kpc, comparable to the stellar 
mass. 

The extension of this MOND model to 1 Mpc predicts a 
circular velocity corresponding to a Newtonian mass of 4.5 x 
10'"^ Mq - which agrees better with the DBA mass. We also cau- 
tion that the MOND dispersion at a radius of about 200 kpc 
has not yet been ruled out (Schuberth et al. .2007) . Moreover, 
the amount of gas beyond 200 kpc is not well-constrained. The 
necessary dark halo does not falsify MOND, but probably re- 
flects the general finding that MOND fails in galaxy clusters 
without invoking a dark matter component. Sanders ( 120071 1 pro- 
poses neutrinos; but since even very cold neutrinos are not ex- 
pected to clump on scales below a few hundred kpc, this would 
not be a viable hypothesis for NGC 1399. Other possibilities are 
molecules (but then how avoid star formation?) or warm-hot in- 
tergalactic gas (e.g. Takei et al. 2007). Given that only 10% of 
the baryons at low redshift have been identified (e.g. Bregman 
120071 ). the amount of gas needed is not in imbalance with the 
missing baryons. It would, however, be surprising to find such 



an amount of gas. Future observations with increased sensitivity 
might well potentially falsify this scenario. 



7. Conclusions 

We have used our new sample of about 660 velocities of GCs to 
constrain the dark halo of NGC 1399 within 80 kpc. The result- 
ing halo, assumed to be of the NFW type, is more concentrated 
than the halos from cosmological simulations, but marginally in 
agreement. A more massive alternative halo, which would ac- 
count for the total mass of the Fornax cluster as estimated by 
Drinkwater et al. (2001 ), may also be consistent with our data. 
We are thus not able to conclusively constrain the cluster-wide 
halo from probes within 80 kpc. We also consider the MOND 
case and find that MOND is not able to reproduce the observed 
GC kinematics under isotropy without invoking an additional 
dark matter component, as already found for other galaxy clus- 
ters. A realistic radial anisotropy does not change this conclu- 
sion. Moreover, only the isotropic model reproduces the X-ray 
mass well. The hypothetical dark halo must have a small core, 
and thus neutrinos do not seem to be viable candidates. It may 
be in gaseous form. Under our assumptions (particularly the in- 
terpolation between the MONDian and the Newtonian regime), 
MOND makes a strong prediction, which perhaps will be falsi- 
fiable in the future. 
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